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' Abstract. In the Wick-Cutkosky model we analyze nonperturbatively, in 

■ light-front dynamics, the contributions of two-body and higher Fock sectors 

! to the total norm and electromagnetic form factor. It turns out that two- and 

three-body sectors always dominate. For maximal value of coupling constant 
a = 2tt, corresponding to zero bound state mass M = 0, they contribute 90% 
to the norm. With decrease of a the two-body contribution increases up to 
100%. The form factor asymptotic is always determined by two-body sector. 



1 Introduction 

In field theory, the state vector \p) is described by an infinite set of the Fock 
' components, corresponding to different numbers of particles. In light-front dy- 

namics 012] the state vector is defined on the light front plane uj-x = 0, where u 
is the null four- vector (to 2 = 0). The wave functions are expressed in terms of the 
variables k±, x: tp = ip(ki±,xi; k2±,X2] ■ ■ ■ ; k n ±,x n ). The total norm (equaled to 
1) is given by sum over all the sectors: ^2 n N n = 1, where n-body contribution 
N n reads: 

N n = (2tt) 3 J \ip(k 1 _ L ,xr,k2±,x 2 ;...;k n ±,x n )\ 2 

In applications, the infinite set of the Fock components is usually truncated 
to a few components only. The belief that a given Fock sector dominates (with 
two or three quarks, for instance) is often based on intuitive expectations and 
on "experimental evidences" rather than on field-theoretical analysis. 

In the Wick-Cutkosky model two massive scalar particles interact by the lad- 
der exchange of massless scalar particles. Two-body sector contains two massive 
particles. Higher sectors contain two massive and 1,2,... massless constituents. 
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In the present paper, based on the work we present the results of our 
study in the Wick-Cutkosky model of contributions of the two- and three-body 
sectors to the total norm. Subtracting them from 1, we get total contribution of 
all the sectors with n > 4. Besides, we also calculate their contributions to the 
electromagnetic form factor. Calculations are carried out nonperturbatively in 
full range of binding energy < B < 2m. 



2 Bethe-Salpeter amplitude in Wick-Cutkosky model 

We use Bethe-Salpeter (BS) amplitude known explicitly in the Wick-Cutkosky 
model [3|- For the ground state with zero angular momentum it reads: 

$( k v) = L_ f +1 9M(z)dz 

{ ,P) V^J-i (m 2 -M 2 /A-k 2 -zp-k-ie) 3 ' { ) 

where k and p are relative and total four-momenta, m is the massive constituent 
mass, M is the mass of the composite system. The representation © is valid 
and exact for the zero-mass exchange. The function gM( z ) 1S determined by the 
integral equation: 

9m(z) = — j ^K(z,z')g M (z')dz' (3) 

with the kernel: 



K(z,z') 



2 

m 

m 2 - i(l - z' 2 )M 2 



Here a = g 2 / (lQirm 2 ) and g is the coupling constant in the interaction Hamilto- 
nian H mt = — gip 2 (x)x{x) ■ In nonrelativistic limit the interaction is reduced to 
the Coulomb potential V(r) = — j. 

The normalization condition for gu^z) is found from the requirement that 
the full electromagnetic form factor Ff u u(Q 2 ) (calculated with full state vector 
\p) and, hence, incorporating all the Fock components) equals to 1 at Q = 0. 
Form factor is expressed in terms of the BS amplitude: 

(p + pTFfulliQ 2 ) (4) 

/<i 4 k ( 1 \ / 1 

{p + p' -2kY (m 2 - k 2 ) <P i^-p - k,pj $ i^-p' - k,p' 

We substitute here the BS amplitude © and find the normalization of gM( z ) 
from the equality Ff u u(0) = 1. The details of calculations are given in [3]. 

The function gM{z) is found from (j3J) analytically in the limiting cases of 
small binding energy B = 2m — M (a^O, B — > 0, M — > 2m) and of extremely 
large binding energy (a = 2ir, B = 2m, M = 0). In the case M — > 2m it reads: 



<IM ( : ) = 8\/27ra 5/2 m 3 ( 1 + — log a 



\z\ + -^(l + \z\)\og(z 2 + a 2 /4) 



(5) 

In contrast to the solution found in [3], eq. © is calculated to the next a order, 
keeping, however, the leading log a term (i.e., neglecting const relative to log a). 
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In the opposite case M = gM=o( z ) has the form: 

g M =o(z) = 6^/30vr 3 / 2 m 3 (l - z 2 ) . 

For arbitrary M the function gM(z) is found from @ numerically. 



(6) 



3 Two- and three-body contributions 

Knowing the BS amplitude, we extract from it the two-body wave function pQ: 



(7) 



This relation is independent of any model. In Wick-Cutkosky model, substituting 
© into (J7J), we find: 



x(l - x)g M {^ ~ 2x) 
2y/n(k 2 , + m 2 - x(l - x)M 2 



2 ■ 



(8) 



Substituting @ into we obtain the two-body contribution to the full 
normalization: 

1 f 1 x{l - x) g 2 M (2x - l)dx 



No 



192 7T 3 Jq 



(m 2 -x(l -x)M 2 y 



(9) 



Three-body contribution N3 is found in [3] by calculating the amplitudes Fig. 
H(at Q = 0), where two-body vertices are determined by the wave function (|SJ). 




Figure 1. Three-body contributions to form factor. 
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4 Results 

For small a, with eq. (J5J) for qm{z), the contributions N2 and JV3 to the total 
norm are found analytically (up to a log a order): 

AT 2 = 1-— logi, N 3 = —log-, N n>i = 0{a 2 ). (10) 
7r a tt a 

For a = 2-7T (B = 2m, M = 0), with qm{z) given by ©, we get: 

iV 2 = JL~ 6 4% ; Ar 3 «26%, iV n > 4 ~10%. (11) 

For a in the interval 2tt > a > 0, corresponding to < M < 2m, the values 
of N2, N3 and A^ n >4 vs. M are found numerically and they are shown in Fig. |2I 




Figure 2. Contributions to the total norm N n =2 + N n =3 + N n >i — 1 of the Fock sectors with 
the constituent numbers n — 2, n — 3 and n > 4 vs. the bound state mass M (in units of m). 



We find that two-body sector always dominates. The sum ./V2+-/V3 contributes 
90% even in the extremely strong coupling case, as we see in (|11|) . This result is 
non-trivial, since for a = 2ir one might expect just the opposite relation of the 
N2 + -/V3 and -/V n >4 contributions. For any a, asymptotic behavior of the form 
factor Ff u n(Q 2 ) is determined by the two-body Fock sector [H]. 
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